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^ ■ Abstract 

Floquet formulation of time-dependent density-functional theory is revisited in light of its recent 
criticism [Maitra and Burke, Chem. Phys. Lett. 359 (2002), 237]. It is shown that Floquet 
theory is well founded and its criticism has overlooked important points of both the Runge-Gross 
formalism and Floquet formulation itself. We substantiate our analysis by examples similar to that 
considered by Maitra and Burke. 
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Density functional theory (DFT) is a well established theory for the ground-states of 
Coulombic systems. It is conceptually simple and practically useful in various branches of 
ohysics, chemistry and material science. Based on the work of Hohenberg and Kohn (HK) 
1] and Kohn and Sham (KS) DFT has become a practical tool for the electronic struc- 
ture calculation of atoms, molecules and solids [3]. Time-independent density_ functional 



formalism is also being extended to deal with excited-states 
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llj. On 



the other hand, time- dependent phenomena are not accessible through traditional DFT. 

is the generalization of the 

ground-state DFT to time-dependent problems. The development of TDDFT is relatively 
recent. The central result of modern TDDFT is a set of time- dependent Kohn-Sham 
equations which are structurally similar to the time- dependent Hartree-Fock equations 
but include in principle exactly all the many-body effects through a local time-dependent 
exchange-correlation potential. TDDFT allows access to lots of interesting and important 
phenomena which can't be easily dealt with in static DFT. It has become popular for 
studying atoms and molecules in intense laser fields 11 51 1 . calculating the excitation spectra 

nn 

and oscillator strengths [16|, H7J, frequency dependent polarizabilities, hyperpolarizabilities 



[181 ] and optical response of solids, etc. 
control problems 



19]. TDDFT is especially very useful for quantum 
20l | . This is because for intense laser fields where correlation effects are 
quite crucial, TDDFT appears to be the only way of studying the quantum dynamics of a 
many-electron system. 



Initial steps towards the rigorous foundation of TDDFT were taken by Deb and Ghosh 
ill ^ and by Bartolotti [ljj. They formulated and explored HK and KS type theorems 
for the time-dependent density. Each of these derivations, however, was restricted to a 
particular class of allowable time- dependent potentials (to potentials periodic in time). 

oped TDDF 



22 



24, 



25| 



Because of the periodicity of the potential in time, Deb and Ghosh deve 
with the Floquet formulation. Further in recent times Telnov and Chu 2l|, 
have developed approaches for the nonperturbative treatment of the strong field processes. 
These are the approaches which involve the transformation of the TDDFT equations into an 



equivalent time-independent generalized Floquet matrix eigenvalue problem. A variational 
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perturbation method based on the Floquet theory has also been developed and applied [26] 
to study optical properties of atoms within TDDFT. This perturbative Floquet approach 
is suitable for the treatment of weak-field processes. A general proof of density to potential 
mapping (p(f,t) — > v(r,t)) for a time-dependent density was given by Runge and Gross 
[jjj]. Runge and Gross (RG) ]14] formally established the TDDFT by showing that for 
a given initial state the evolving density p(f, t) uniquely determines the corresponding 
time-dependent potential v(r,t). Two features of RG formalism are noteworthy: one that 
the time-dependent wavefunction for a given density has functional dependence on both the 
density as well as the initial wavefunction \& • Secondly the formalism is valid for potentials 
that are Taylor series expandable at the initial time t = to- This makes the formalism appli- 
cable to suddenly switched-on potentials but not to the adiabatically switched-on potentials. 

Like in the case of static DFT, in TDDFT also an interacting many-electron system is 
mapped onto a fictitious non-interacting Kohn-Sham system with the same time- dependent 



density as the interacting system. In Floquet formulation 12J, [l3[] the Kohn-Sham system is 
developed in terms of the periodic time- dependent density only. On the other hand, in the 
the RG formalism the Kohn-Sham wavefunction is a functional of both the time-dependent 



density and an initial state 27|, |28[]. The uniqueness of the KS system follows from the 
one-to-one p{r,t) — > v(r,t) mapping. However, in a recent work [29| by Maitra and Burke 
it has been suggested that a KS system within the Floquet formalism is not unique. The 
non- uniqueness of the KS system has been attributed to its functional dependence on the 



121 ] may not 



initial state. Further, they suggest that the Hohenberg-Kohn like theorem 
exist in Floquet theory of TDDFT because the "ground Floquet state" cannot be defined 
properly. This appears to have created some doubts 3(3] about the Floquet theory in 
TDDFT. It is the purpose of this paper to analyze the Floquet theory of TDDFT and show 
that it is perfectly valid. The present analysis also brings out subtle differences between the 
RG and Floquet theory and shows the two formulations to be distinct and founded on two 
different principles. As such Floquet theory of TDDFT cannot be thought of as a particular 
case of the RG theory and the RG arguments can't be applied to Floquet formulation. We 
begin with a short description of the foundations of Floquet and RG theory in TDDFT. 
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For a time-dependent (TD) Hamiltonian H(t) periodic in time, i.e. 



H(t) = H(t + r) , 



(1) 



where the period r is positive with uj = — , the TD Scrodinger equation (atomic units are 
used throughout) 



has solutions 



31 



#„(f,t) = o , 



32] of the form 



* n (r, t) = n (f, t) e ie ™*; n (f, t) = <f> n (r, t + r) , 



(2) 



(3) 



where the time-periodic functions <f) n (r,t) are termed as the quasienergy eigen states (QES) 
and e n are referred to as the the quasienergy eigenvalues. The steady-state wavefunctions 
4> n (f,t) are solutions of 

d 



H(t) - i 



dt 



4> n (r,t) = e n (j) n {r,t) 



(4) 



where n (r,t) are square integrable and e„ are real numbers. The quasienergies, 
Eq, £i, .., £j, ... are defined 'modulo (a;)'. However, they can be ordered in increasing or- 
der Eo,£i,... by making sure that as the strength of the applied time-dependent potential 
goes to zero, each quasi-energy go to its unperturbed counterpart [33j. Thus as long as to is 
not one of the resonant frequencies, there exist the "ground-state energy" e and "excited- 
state energies" ei, e n in the steady-state formalism (see Appendix A for details). The 
theory of solutions in time-periodic Hamiltonian can be thought of like the stationary-state 
theory but in an extended Hilbert space (R + T) that includes in addition to the space 
dependent function, time-periodic functions also. The operator 



H(t) = H(t) - , 



(5) 



is called the Hamiltonian for the steady states in the composite Hilbert space, which re- 
sembles in many way the Hamiltonian for bound-states. The scalar product in this space is 
defined as 
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{(</#>} = \J dtj f{r,t)ijj{r,t)dr (6) 

i.e. in addition to the space integral, an integral over time is also taken. Here the curly 
brackets '{ }' indicate the time average over a period. Thus the quasienergy functional is 
given by 

e n [0 n ] = {(0 n |W(i)|0 n >} • (7) 

We point out - and this is important from TDDFT point of view - that the steady-state solu- 
tions are obtained by an adiabatic switching of the periodic pot ential. Finally applying the 



variational principle, each state (f> n (r, t) can also be obtained 32] by making the expectation 
value |^0 n H(t) — i-j^ 0n^j stationary. The stationarity principle becomes a minimum 
principle for the "ground-state" (Here the "ground-state" refers to a steady-state having the 
lowest quasienergy.) by assuming completeness of the set {<fi„} over the (R + T) space. The 
general proof of the minimum principle for the ground-Floquet state is given in Appendix A. 

Using the minimum property of the ground-state quasienergy functional a Hohenberg- 
Kohn like theorem can be proved [3] for time-dependent periodic densities also. We repro- 
duce the proof below. Let us consider the nondegenerate ground-state of a system charac- 
terized by the Hamiltonian H(t). Suppose there exist two different Floquet ground-states 
<f) and <p' corresponding to the steady state Hamiltonians 7i and TC', or more specifically 
corresponding to the external potentials v(f,t) and v'(f,t), both of which give rise to the 
same TD density p(r, t). Using the minimal property for the ground-state we get 



1 



t2 



e' < e + - J dt j [v'(f, t) - v(f, t)}p(f, t)d s r . (8) 
Similarly for the ground-state of the Hamiltonian TC 



1 



T 



*2 



tl 



£<s' + - dt [v(f, t) - v'(f, t)}p(f, t)d 3 r . (9) 



Adding Eqs. (jHJ) and (El) results in the absurdity, 

e + e' <e' + e . (10) 



Hence two different potentials v(f,t) and v'(r,t) cannot give the same density p(r,t), 
which implies that the time-dependent potential v(f,t) is a unique functional of the 
time-dependent density p{f,t). Accordingly, TC, H and and indeed any "ground-state" 
property are all unique functional of p(f, t) . Thus for Floquet states the HK theorem is [12|] 
: "the density corresponding to the ground-state quasienergy of a time-periodic Hamiltonian 
determines the corresponding external time-periodic potential v ext (f, t) uniquel y" . For 



example in perturbation theory up to the second order in energy, is minimum [26( with 
respect to p^ for frequencies less than the first excitation frequency (i.e. uj < uj 10 = Ei—E , 
where E and E x are the unperturbed ground and excited-state energies). Thus in the 
linear response regime, the theorem would apply to p^ up to freque ncy uj < com- This 



theorem is then the foundation of Floquet theory of TDDFT [2l|, [22, l23|, |24|, l25|, |26|]. We 
point out that the theory does not have any initial state dependence but requires adiabatic 
switching of the applied external potential. Assuming the v— representability, an equivalent 
Kohn-Sham system for an N— electron system can also be developed. 

As pointed out earlier, in RG theory [ijj] the proof of the density-to-potential mapping 
for time-dependent systems is based directly on the TD Schrodinger equation. In this 
theory , it is assumed that the time-dependent potential v (r, t) is turned on at a particular 
time t = to and all systems are taken to evolve from the same initial state wavefunction 
^(to) = ^o- It is then shown that if the potential has a Taylor series expansion around 
the initial time t = t , then the time-dependent density p(r,t) determines the potential 
v(f,t) uniquely. In turn ty(r,t) is also determined by p(r,t). Thus the RG theory is 
quite distinct from the Hohenberg-Kohn theorem or its steady-state DFT counterpart 
(Deb-Ghosh theorem proved above), which are based on the minimization of energy (or 
quasienergy in the Floquet formalism). Notice that because of the Taylor series expansion 
requirement of v(f, t), RG theorem is applicable to suddenly switched-on potentials but not 
to adiabatically switched-on potentials. Table I gives a comparison of the two theories. 

It is clear from the comparison given in Table I that the two theories are quite distinct 
and their domain of application is also different. As such Floquet formulation of TDDFT 
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TABLE I: Comparison of Floquet theory and RG theory of TDDFT 



Floquet Theory 



• Based on the minimum principle 
of the ground-state energy. 

• Steady-state is reached by 
adiabatic switching of the 
time-periodic potential. 

• v(f,t) is functional of p(r,t) only. 



RG Theory 



• Based on the evolution of 
wavefunction from t = to. 

• Wavefunctions are obtained 
from potentials analytic a t = to- 
Adiabatic switching is ruled out. 

• v(f,t) is functional of p(r,t) as 
well as ^o- 
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13 



21 



22 



23 



24 



25 



261 ] is not a particular case of the RG theory. In light of this 



291 ] and then go on to comment on 



we now analyze the paper by Maitra and Burke 
what does the distinction between the Floquet and the RG theory mean in the context of 
calculation of quantities like polarizabilities a(u) and excitation energies etc. using TDDFT. 



Maitra and Burke consider the example of two non-interacting electrons in a one- 
dimensional harmonic oscillator potential subject to an adiabatically switched-on time- 
dependent potential periodic in time. Thus the Hamiltonian they consider is 

Id 2 1 

H = — - + -ujIx 2 + \xSin(ut) . (11) 

2 ax 2 2 

Due to the lack of many exactly known analytic solutions of time-dependent systems, we 
also consider the same system. The QESs for the above mentioned Hamiltonian are known 

n 

analytically [34J to be 

cj) n (x,t) = i; n (x(t))e {i ^ AxCos ^ t)+a{t)} , (12) 

where x/j n are the eigenstates of the static Harmonic oscillator and x(t) = x — ASin(ut). The 
parameters of the wavefunction are given as 



and 

f Sin(2uit) 



«(*) = A 21 — — ^ 2 , • ' (14) 

The corresponding quasienergy for the QES 4> n (x, t) is 

(n + -)ujq + 



+ Cosioot) - 1 - u 2 Cos(ut) x 



moduloiuj) (15) 

We first point out that the quasienergies have been expressed in such a manner that as 
A — > 0, the energy eigenvalues go to their respective time-independent eigenvalues. Thus 
energy defined ' modulo (u;)' does not create any difficulty Secondly, as long as u ^ u , there 
is clearly a well defined "ground-state" in these solutions. Thus so = ^ + 77 T t\ represents 
the "ground-state" energy of the system, whereas the other energies are the "excited-state" 
energies. 

Maitra and Burke consider a singlet state with one electron occupying the ground-state 
(i.e. n — 0) and the other the first excited-state (i.e. n = 1) quasienergy orbitals. This 
gives the density of the system to be 

p(x, t) = ^ (1 + 2u x(t) 2 ) e-^ 2 (16) 

Now they generate the same density by another system with a different periodic potential. 
For doing this they consider a Floquet state which is also a spin singlet but with one doubly 
occupied steady state. Thus 

$(x u x 2 ,t) = 0(x u t)0(x 2 ,t)e 2m , (17) 

with 

0(x,t) = y^e*<*'*> , (18) 

where (3(x,t) is a real time-periodic function, /3(x,t + r) = f3(x,t). By inverting the TD 
Schrodinger equation one will have a different potential 

, . 10"(£,t) <P(x,t) . . 

v(x,t) = - ~ ' +i-^- L +e 19 
2 0(x,t) J>(x,t) 



8 



In this s is the orbital quasienergy for the second Floquet state. The function (3(x,t) is 
determined by assuming the potential and quasienergy to be real. Now the new potential 
and phase for the second Floquet state are given by 

s(m) = ^im 2 {1 - (1+2 4 wy - rra^y } + 

ujq u 2 A 2 Sin 2 (ut) 
1 + 2uj x{t) 2 +u ° 2 

(3(x, t) = AuCos(cot)x -A 2 (lu 2 -4) Sm( : 2Ujt) (20) 

V 2 J 4:U> 

The phase fi(x,t) is uniquely determined up to a purely time- dependent function. At 

2 2 

large x both the potentials v(x,t),v'(x,t) — > — h XxSin(cot) and the quasienergy is 

Based on the fact that they have been able to generate one density from two different 
potentials, Maitra and Burke then conclude that this is a manifestation of the initial 
state-dependence of TDDFT. Thus they find Floquet theory of TDDFT to be incomplete. 
We now show that Maitra and Burke are incorrect on two counts. First, they have 
considered a system in steady-state solution that is supposed to have been obtained by 
an adiabatic switching. Thus RG theory is not applicable to this system. Therefore 
one cannot talk about the dependence on the initial state which is a specific feature 
of the RG theory. The problem falls in the domain of steady-state solutions and as 
such should be looked at within the Floquet theory of TDDFT. As discussed earlier, 
the domain of applicability of the Floquet theory of TDDFT is precisely the kind of 
example that has been considered above. Secondly, the Floquet theory of TDDFT is 
similar to the stationary-state DFT and as the derivation earlier shows, it is applicable 
to only the ground-states of the steady-states. Maitra and Burke on the other hand have 
applied it to an excited-state of the steady-state solution and shown that an excited-state 
density can be generated by two different potentials. That however, does not invalidate Flo- 
quet theory of TDDFT. This point needs further elaboration and we do that in the following. 

For our discussion, we show that even in the stationary-state theory, we can generate 
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FIG. 1: Shown in the figure are the two potentials (lower panel) giving rise to the same excited- 
state density (upper panel)of a harmonic ID oscillator by putting one electron in the ground and 
other electron in the first excited-state. 

a given excited-state density by two different potentials , Again consider 

two electrons in a one-dimensional harmonic oscillator potential. If the two electrons are in 
the ground-state, the potential is determined uniquely by the Hohenberg-Kohn theorem. 
However, for excited-states, there is no such theorem and more than one potentials can give 
the same density. For this let us consider one electron in n = and one in n = 1 state. The 
corresponding density is shown in the upper panel of Fig. 1. One potential corresponding 
to this density is obviously the harmonic potential v(x) = itself. Another potential 

can be generated by putting both the electrons in the ground-state orbital and the n by 



inverting the Schrodinger equation or by using the van-Leeuwen-Baerends method 



35) 
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FIG. 2: Two potentials (lower panel) yielding the same excited-state density (upper panel) for an 
excited state of the two-electron in an ID harmonic oscillator. 



These two potentials are shown in the lower panel of Fig. 1 . We consider one more example 
similar to that discussed above, where the excited-state density of two non-interacting 
electrons moving in a one-dimensional harmonic oscillator potential is obtained by putting 
one electron in n = and the other in n = 2 state. Now this density (shown in the upper 
panel of Fig. 2) is produced by an alternative potential, which is shown in the lower panel 
of Fig. 2 along with the original potential. Thus if a system is in excited-state, more 
than one potentials can give the same density. On the other hand for the ground-state, 
occupation is uniquely defined and so is the potential. More examples of this kind exist in 



the literature 



a, y, y, 



10 



11] . Exactly in the same manner as in stationary-state theory, 



in Floquet theory of TDDFT also the "ground-state" gives the potential uniquely but for 
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"excited-states" more than one potential can be found. This is precisely what has been 
done in (29]. Maitra and Burke do state in their work that their results resemble those of 
excited-states in time-independent problems but fail to make further connection with the 
latter. 



Having explained the work of Maitra and Burke, we now also comment on the TDDFT 
calculations in light of the Floquet and RG theories of TDDFT. It is clear that in calculating 
frequency dependent polarizability a(u>), it is the steady-state of a system that is employed. 



Thus in ca 
strengths, 



cu 



ation of at(u), and related quantities such as excitation energies and oscillator 
16l ] it is the Floquet theory of TDDFT rather than the RG theory that is being 
applied. Finally, one may raise a question if Floquet theory of TDDFT is applicable only to 
"ground-states". The answer is that even for steady-state "excited-states", a theory similar 



to the stationary-state excited-state DFT 
is not our main concern here. 



tor steaay-st, 

ifl Sells, 
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ll| can be developed but that 



To conclude, we have shown that Floquet theory of TDDFT is well founded and is 
distinct from the RG theory. Further, its recent criticism by Maitra and Burke is easily 
explained on the basis of a careful analysis of Floquet theory and RG theory. 
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APPENDIX A 



In this appendix, we show that if one confines the quasienergies to a range such that 
they go to their unperturbed counterpart when the time-periodic potential is turned off, 
then there is a well defined ground-state quasienergy satisfying the minimum variational 
principle. We first point out that defining the eigenenergies as suggested above is equivalent 
to taking the corresponding quasienergy state such that it has no free time-dependent factor 
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of the form \ J2 p e ipu)t j , where {p = 0, ±1, ±2, ±i, ..}, associated with it as the strength 
of the time-periodic potential goes to zero . Let us denote such quasienergy states as 
<f>o(r, t), 0i (r, t), 0i(f, t), ... with the corresponding eigenvalues £o, £i, £j, ... such that 
£ < £i < ... < £, < .... The quasienergy states satisfy the normalization condition 

/ <pl n (f,t)tp n (f,t)df= S mn , (Al) 

as is easily seen by the hermiticity of the Hamiltonian coupled with its adiabatic switching. 

Let us consider a Hamiltonian H(f,t) = H (f) + v(f,t) with v(f,t) being time periodic. 
Now a trial periodic function $(r, t) can be expanded in terms of unperturbed states 0i(r)' s 
of #o(r) and its associated states 4>i(r) e rpuJt as 



$(r,t) = a 4> (r)+a 1 4> 1 (r) + .... + a i 4> i (r) + 




In the expansion above, each coefficient a\ puJ> (p ^ 0) should become zero as v(f,t) — > 0. 
Otherwise coefficient af^ would have two components: one arising from the applied time- 
dependent potential causing unperturbed states to mix and the other is the coefficient of 
physically equivalent state (f) i e' ipuJt . To see this let us add $(r, t)e luJt to $(r, t) to get 



$'(r, t) = a (po(r) + a X (j)i(r) + .... + a^r) + 

+ (4 W) + «o) 0o(r)e^ + (a< w) + a ) 0i(r)e^ + .... + (a, H + a ) ^e*" + 

= (l + e^)$(r,t) (A3) 

For normalized wavefunctions the factor in front of $(r, t) gives rise to a purely time- 
dependent phase factor and can therefore be ignored. This is best illustrated if we look 
at the unperturbed time-independent problem in the Floquet formulation. 
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In the case of time-independent problem in Floquet formulation the time-dependent wave- 
function corresponding to a state <f>i(r) can be written in two equivalent forms as 



-iEit 



-i(E+pw)t 



(A4) 



where 4>i(r,t) are the Floquet states satisfies 







H(t) - i— ) Mr, t) = (Ei + pufyiir, t) . 



(A5) 



However, each (j>i{r, t) gives the same ^(r, t). So they are all physically equivalent. If we 
were to expand a trial Floquet function in terms of {4>i(r}} , {cj)i(r)e luJt } , {<j)i(r)e ipu}t }, all 
it will do is give a complicated time-dependent phase factor in front of (f>i(r). How does this 
physical equivalence gets reflected variationally is discussed next. 



Let us take a trial periodic wavefunction as 



3 II ] for the unperturbed problem as 



$(r,t) = ^a l (p V l (r1e 



ipujt 



ip 



with all 0's and their equivalent states included. The quasienergy of the system is 



1 

T 



$(f,t) 



H(t) - i 



d_ 
di 



$(f,t) ) dt 



\ f e-^Ej + q^af faftt) e^dfdt 

^SijSpqiEj + qu)a\ 

^\a?\E i+ pu) . 



(p)* (<?) 



Since $(r, t) is normalized at all times, to obtain af \ the functional above is made stationary 



(A6) 



(A7) 



261 ] with the condition 



i I ($(f,t)|$(f,t)) = E 

ip 



a] 



(p) 



1 . 



(A8) 



Using the techniques of Lagrange multipliers this leads to the equation 



a^ ) [E i +pu-fi] = , 



(A9) 
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where p, is the Lagrange multiplier and gives the quasienergy of the system. The equation 
above gives p = Ei + pu with of 7^ but all other of s vanishes. Thus stationary 
variational procedure picks Floquet states from only one particular zone {i.e. 'p' is fixed) to 
represent the system. However, irrespective of which 'p' we take, the final wavefunction is 
the same. Thus Floquet states from different zones do not represent different states. As such 
one may restrict oneself to only one zone in the expansion of the Floquet state. Question 
that arise is: does restricting oneself to one particular zone lead to a minimum principle? 
Actually, it does as we will now show. It is best shown for p = zone but the result is true 
for any 'p'. Thus 



$(r,t) = $>i°Vi(*0 , (A10) 

i,p=0 



d 



<*i£-^i*> = EH 0) £ * > EK 0) £ o , (ah) 

i i 

by the standard variational argument. 

Having discussed the time-independent case, we now discuss the time-dependent case. In 
this case the basis functions are the Floquet state (pi(f,t) of the full Hamiltonian and their 
physically equivalent counterparts <pi(f,t)e ipu}t . If we expand a trial Floquet state $(r, t) as 

<S>(r,t) = Y,c ( fUi(r,t)e^ t , (A12) 

ip 

we again argue that including Floquet states from all zones does not really give us any new 
information. So the expansion should be restricted to only one zone. Mathematically it is 
again shown as was done above. The approximate quasienergy is given as 

f J(®\H-*§- t \m = J2\ c ? ) \ 2 ^ + P^ ■ (A13) 

n 

So the stationarity [26( of the time averaged expectation value in Eq JA13l under the 
constraint ~ J dt = 1 gives 

cf ] {ei + pu-E} = (A14) 
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implying again that E = £j + pu with cf' - only in one zone. However, Floquet states 
from each zone give the same wavefunction ^ = e _ ^ £+pa, ^0i(r, t)e ipu)t . Thus in expanding 
an approximate $(r, t) one can restrict oneself to one particular zone. This leads to a 
minimum principle following the standard arguments. 

The question arises how do we make sure that the trial wavefunction comprises Floquet 
states from one particular zone only. This is best done for p = zone by making sure that if 
we take v(r,t) = 0, the trial wavefunction must be time-independent. If Floquet states from 
other zones are also present in the trial wavefunction, it will not become time-independent as 
the time-dependent potential becomes zero. With such a trial wavefunction the quasienergy 
follow a minimum principle as shown above. 
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